Shot noise and Coulomb effects on non-local electron transport in 
normal-superconducting-normal heterostructures 
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We argue that Coulomb interaction can strongly influence non-local electron transport in normal- 
superconducting-normal structures and emphasize direct relation between Coulomb efTects and non- 
local shot noise. In the tunneling limit non-local differential conductance is found to have an S-like 
shape and can turn negative at non-zero bias. At high transmissions crossed Andreev reflection 
yields positive noise cross-correlations and Coulomb anti-blockade of non-local electron transport. 



I. INTRODUCTION 

Discreteness of electron charge has a number of funda- 
mental physical consequences, such as, e.g., shot noise in 
mesoscopic conductors^ and Coulomb blockade of charge 
transfer in tunnel junctions^. About 10 years ago it was 
realized'^'* that these two seemingly different phenom- 
ena are closely related to each other: Coulomb blockade 
turns out to be stronger in conductors with bigger shot 
noise. This fundamental relation was subsequently con- 
firmed in experiments^. A close link between shot noise 
and Coulomb blockade exists not only in normal con- 
ductors but also in hybrid normal-superconducting (NS) 
structures^, where doubling of elementary charge due to 
Andreev reflection becomes important. 

Can the above relation be further extended to in- 
clude non-local effects? A non-trivial example is pro- 
vided by normal-superconducting-normal (NSN) sys- 
tems where entanglement between electrons in different 
normal terminals can be realized. Non-local electron 
transport in such systems is determined by an interplay 
between elastic cotunneling (EC) and crossed Andreev 
reflection (CAR) and was recently investigated both 
experimentally^"— and theoreticallyi^rJi (see also fur- 
ther refs. therein). While non- interacting theory predicts 
that CAR never dominates over direct electron transfer 
(hence, no sign change of non-local signal could occur), 
both positive and negative non-local signals have been 
detectcd^.i^ d°di . It was argued that CAR could prevail 
over EC in the presence of Coulomb interactionsi^ or an 
external ac field-^. Negative non-local conductance was 
also predicted in interacting single-level quantum dots 
in-between normal and superconducting terminalsii. 

Despite these developments no general theory describ- 
ing the effect of electron-electron interactions on non- 
local transport in NSN structures was available until 
now. Below we will construct such a theory and demon- 
strate that interaction effects in non-local transport and 
non-local shot noise in such systems are intimately re- 
lated. This relation, however, turns out to be much more 
subtle than in the local case^i^i^ merely because of (a) 
a variety of different processes contributing to non-local 
shot noise and (6) positive cross-correlations which may 
occur in normal-superconducting hybrids^i^ (in contrast 
to normal conductors where cross-correlations of fluctu- 
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FIG. 1: Schematics of the system under consideration. 



ating currents are negative^). In tunnel NSN systems 
EC and CAR provide respectively negative and positive 
contributions to non-local shot noise^^-^". Here we will 
analyze non-local shot noise beyond the tunneling limit 
and find that at higher transmissions also direct elec- 
tron transfer can yield positive cross-correlations in ad- 
dition to CAR. At full transmissions only positive cross- 
correlations due to CAR survive and yield Coulomb anti- 
blockade of non-local electron transport. 

The paper is organized as follows. In Sec. 2 we describe 
our model and derive an effective action for NSN system 
under consideration. In Sec. 3 we formulate the Langevin 
equations describing real time dynamics of fluctuating 
voltages and currents and derive the general expressions 
for both local and non-local current-current correlators 
describing shot noise in our system at arbitrary barrier 
transmissions and arbitrary frequencies. Sec. 4 is de- 
voted to the effects of Coulomb interaction on both local 
and non-local conductances of our NSN device. A brief 
summary of our key observations is presented in Sec. 5. 
Some technical details are outlined in Appendix. 



II. THE MODEL AND EFFECTIVE ACTION 

We will consider a hybrid structure consisting of two 
normal electrodes coupled to a superconductor via NS 
barriers with local subgap conductances Gn and G22 and 
capacitances Ci and C2 (Fig. 1). External voltages Vi 
and V2 are applied to normal electrodes via Ohmic shunts 
with conductances Gf^ and G^. Weak electromagnetic 



coupling between two NS barriers (e.g. via modes propa- 
gating in the superconductor—) will be disregarded. The 
Hamiltonian of the system reads 



H — Hi + H2 + Hs + Hta + Ht,2i 



(1) 



where 



Hr= J2 '^^'^la 



a=t,i 



V2 



1,2, 



are the Hamiltonians of the normal metals, m is electron 
mass, fj, is the chemical potential, 



Hs 



dx 



E^"'^ 



-^^ ("2^ " '") Xa + Axjxl + A*XiXt 



is the Hamiltonian of the superconductor with order pa- 
rameter A and 

HT,r=Ar Yl Vre'^-i^lxr,o.+Ke-'^-xlJp] (2) 

are tunneling Hamiltonians describing transfer of elec- 
trons across the contacts with area Ar and tunneling 
amplitude tr- For the sake of simplicity we will as- 
sume that both N S barriers are uniform implying that all 
Nr — kpAr/4:Tr conducting channels in the r-th barrier 
are characterized by equal transmission values 



Tr = 4n^iyriys\tr\^/i'^ + TT^iyr'^s\tr\^) 



(3) 



where lyj {j = 1,2, S) is the density of states in the corre- 
sponding terminal. Accordingly, in the low energy limit 
to be considered below local subgap conductances are de- 
fined aS^'^ Grr = {2e^/TT)NrTr, whcrC Tr = T^ / {2 - Trf 

are effective Andreev transmissions of N S barriers. Fi- 
nally, we note that fiuctuating phases Lpr introduced in 
Eq. ^ are linked to the voltage drops across the barriers 
Vr by means of the standard relation ipr — evr and are 
treated as quantum operators. 

As usually, we eliminate fermionic variables and ex- 
press the kernel J of the Keldysh evolution operator via 
path integral over the phase fields^ 



J= f Y[ V^r'^^r(^MiSenA^]+iSTy]), 

<^ ^ 1 n 



(4) 



r=l,2 



where ipr and ipr a-re fiuctuating phases defined re- 
spectively on the forward and backward branches of the 
Keldysh contour, 5cnv is the action of electromagnetic en- 
vironment and the term iSx accounts for electron trans- 
fer between the terminals. In the case of linear Ohmic 
environment considered here one haa^ 
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(eVr-iPrK-Cr'Pr+GfiPr) 
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(5) 



where ipr = {^r + '^f )/2, Vr 
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The term iSt reads 
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(6) 



where 4x4 matrices G represent the inverse Keldysh 
Green functions of (isolated) normal {j = 1,2) and su- 
perconducting (j — S) terminals and tr is diagonal 4x4 
matrix in the Nambu - Keldysh space 
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After some exact manipulations we obtain 
iSr = tr In [l - f\GihGs - 4^2^265] 



(8) 



While the expression ([8]) for the action remains formally 
exact it is still too complicated to be directly employed in 
our calculations. In order to proceed we will make several 
additional steps which yield necessary simplifications. 

As a first step, we restrict ourselves to the limit of high 
conductances 



5r 



2^(Gf + G„)/e">l, 



(9) 



in which case phase fluctuations are weak and it suffices 
to expand the action ([8]) to the second order in (p~ , cf., 
e.g.)^i2i. Technically, we first expand the matrices U to 
the second order in the quantum phases (p~, i.e. we make 
a replacement 



ir^iriipr)[l-l^A-^i 



(10) 



Here irifr) is defined by Eq. ^ with LPr'^ being re- 
placed by the classical phase tpr, and A is a diagonal 
matrix with non-zero elements An = 1, A22 ~ —1, 
A33 = — 1 and A44 — 1. Accordingly, we can write the 
product i\Grir in the form 

ilGrir ^tr + S^.^) + S^' + 0{ip-^), (11) 

where we defined the self-energies 

tr = il{^r)Grir{Vr), (12) 

-[pr'A^tr], 



sp) 



S(2) = --^{ip^^l,tr] + \iPr'A^rVr'A- 



In order to evaluate S^ we employ the Keldysh Green 
functions of the normal leads 

Gr = (G^ + G^) ® d,/2 + (Gf - G^) ® Q.a./2. (13) 

Here the 2x2 matrices G^'^ and Qr are defined as 

nRA 



G^'^ 
G«^^+ 



1 — 2nr 2nr 

2 — 2nr 2nr — 1 



where n^ is the quasiparticle distribution function in the 
r-th normal lead. In equilibrium it coincides with the 
Fermi function n^ — 1/(1 +exp(_E/r)). Neglecting the 
proximity effect in the normal leads and performing the 
summation over the corresponding electron states, we ex- 
press the zeroth order self-energies (J12p in the form 










(14) 



with r^ = 2Tn'r\tr\'^. The function hr{r) in this expres- 
sion differs from zero only at the interface of the r-th 
junction and it obeys the following normalization condi- 
tion 



(frhr{r) = Nr. 



(15) 



For the sake of simplicity, in what follows we will as- 
sume that the barrier cross-sections remain sufficiently 
small and put hr{r) = NrS{r — r^). This assumption 
just implies that each of the barriers has Nr conducting 
channels with identical transmissions Tr ([3J as we already 
indicated above. In this case we can reduce the full coor- 
dinate dependence of the Green functions to that on the 
two indices i and j which label the barriers and, hence, 
can take only two values 1 and 2. Accordingly, e.g., the 
Green function Gs{r,r') reduces to the 2x2 matrix in 
the "junction space" Gg . In addition we should bear in 
mind that G^ are the matrices in the space of conducting 
channels. 

We note that the above assumptions are not really re- 
strictive since they do not affect the general structure of 
our effective action to be derived below. At the same 
time they allow to establish relatively simple expressions 
for the parameters entering in the action. Expanding 
the action (|S]) in powers of ip^ we arrive at the following 
expression 
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where we define the operator 

k= [l-EiGs-EzGs]^ 



(16) 



(17) 



Our second step allows to establish an explicit expres- 
sion for the operator p7)) . Namely, in the interesting for 
us low energy limit 



T,eVr<S: |A| 



(18) 



we can set the energy argument E in the superconduc- 
tor Green function Gs{E) equal to zero. After that Gs 
reduces to the time/energy independent matrix 



Gs{E = 0) = {G§ + Gj) ® {a,/2) 



(19) 



where we introduced the retarded and advanced Green 
functions of the superconductor Gg and Gg . In the limit 
_E ^- they are equal to each other both being 4x4 
matrices in the Nanibu "junction" space 
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(20) 



In addition, each of the matrix elements in Eq. (I^Hl) is 
itself a matrix in the channel space. For instance, G{\ 
and G^2 ^^^ respectively Ni x A^i and A^i x 7V2 matrices, 
while matrices G|\ and G22 have dimensions N2 x A^i 
and N2 X 7V2 respectively. 

As the Keldysh Green function P^ depends neither 
on time nor on the quasiparticle distribution function, 
it commutes with the phase factors e^"^'''-*-' entering the 
self-energies. This observation combined with the multi- 
plication rule for the Q— matrices, Q1Q2 = 1 — Qi + Q2, 
allows us to express the operator K (fTT]) as a linear com- 
bination of these matrices: 



K{t, t') - [(Kr + Ka)/2] ®\5{t- t') 



~i{T^/2)KrXiKa^ 
-i{T^/2)KrX2Ka^ 
-i{T2/2)KrX^Ka^ 

-i{T2/2)KRXiKA^ 



iipi{t) 



-itpiit) 



c^zQi{t-t')aze 



(21) 
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\ipi{t') 



JV2{t)~ 



'U.p-«V2(t') 



'iip2{t) 



(^zQ2{t-t )crze 
c^zQ2{t~t')aze 



-iV2{t') 



Here 4x4 matrices X^ are defined as follows 
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while 4x4 matrices Kr^a read 


Kr = Al = 


1 + i^(Xi + X2) + iy (X3 + X4) 



(22) 



Having established the expression for the operator K 
we can now substitute it into the action ([T6|) . Proceeding 
along these lines and going trough a straightforward but 
rather tedious calculation we arrive at the result which 
still turns out to be too complicated for our present pur- 
poses. Further simplification amounts to neglecting local 
interference terms involving the products of the Green 
functions Gf^, F/J, G22 and F^. Technically this step 
is equivalent to replacing these Green functions by their 
averaged-over-disorder values, which read 



(Gfi) = {Gg) = 0, (F«) = (Fi) 



-TTiys- 



(23) 



After this step we immediately arrive at the central result 
of this section 



iSr = iSii + iS22 + 1812, 



(24) 



where 



iSu 



~^~T I dttpitf^ 



2e^ 



iSi2 = i—^ / dt{tfiLp2 +<P2<Pi 



dtdt 



,Vl{t)S'^2V2{t') 



(26) 



and the term 1^22 is obtained by interchanging the indices 
1 f^ 2 in Eq. (gS]). The functions 5**' read 

S\{ = GiiAf(t-i')(l-A+/^icos[2(^**']) 
+ 2Gi2M(t - t')(ai - '71 cos[2^**']) 
+ (Gi2/2)Af(i- i')(4cos[</7f + ./jf ] 
+ Kj;cos[(^f -vj"']), 
-Gi2Af (t - t')(l - A + A cos[2(pf' ]) 
-Gi2M(t - t') (1 - /32 + h cos[2v?**']) 
+ (Gi2/2)M(t - i')(7+ cos[(^**' + (^**'] 



(27) 



ctt' 
'-'12 



-7_cos[(^f -¥'"]) 



(28) 



Here we denoted (p** — fr{t) — frit') and introduced 
Andreev Fano factors (3r — 1 — Tr- 

The zero bias non-local subgap conductance G12 as 
well as the parameters ar^i^r, Kr,"f are expressed as traces 
of certain combinations of the matrices Kr, Ka and Xj, 



as described in Appendix. In order to reduce them to a 
tractable form we further assume that normal state resis- 
tance R^ of the superconducting wire segment of length 
equal to the superconducting coherence length ^ remains 
much smaller than normal resistances of the barriersi^, 



I.e. 



e'^NrTrR^/ir < 1. 



(29) 



This condition is usually well satisfied for generic sys- 
tems. Eq. (1^^ enables one to treat the Green func- 
tions connecting the two junctions as small parameters 
and expand the traces in Eqs. (jAll) - (jA3l) in powers of 
Gf2 , G21 , F^ and F^ . Keeping the leading corrections 
ex Gf2G|\ and ex F^F^ and making use of the fact that 
G12G21 becomes equal to F(^F^ after averaging over dis- 
order, we arrive at explicit expressions for the parameters 



rjr = 

7± = 



Tr{l - 2Tr)/^TlT2, 
2Tr{l - Tr)/^TiT2, 

±(4t, - 3) + 1/Vn^ ('• = 1,2), 

±1 + (1 - 2ti - 2t2 + 4TiT2)/V^r^, (30) 



and for zero bias non-local conductance 



Gi 



2e^NiN2TiT2 



r<R r<R 

'-'12'-^21- 



(31) 



In the case of diffusive superconductors one has to fur- 
ther average Eq. (j31l) over disorder. In a simple quasi- 
one-dimensional geometry this averaging yields^^ G12 = 
GiiG22Ri e^'^/^ /2, where d is the distance between two 
NS barriers. 

It is important to emphasize that all order terms in t^ 
are fully accounted for in Eqs. ([25|) -([30 |) . i.e. our action 
applies for arbitrary transmission values Ti^2 (or T1.2) 
ranging from zero to one and - similarly to N S systems^ 
- entering in the combination T^ / (2 — Tj.)"^ representing 
Andreev transmissions Tr- 



ill. LANGEVIN EQUATIONS AND SHOT 
NOISE 



The quadratic in ip-^ 2 action is equivalent to the 
Langevin equations^ir— 

Civi + (Gf + Gn)«i - Gi2t'2 = GfVi + £,t + Ci, 

G2^2 + {Gf + G22)«2 - Gi2t'l - GfV2 + ^f + 6,(32) 

which describe the current balance in our system. Here 
^f^ are stochastic variables with pair correlators 



{if{t)if{t'))^GfM{t-t') 



(33) 



describing Gaussian current noise in the shunt resistors, 
while the variables ^ with the correlators 



{umit')) 






(34) 



describe shot noise in NS barriers. Note that we ignore 
the electro-magnetic coupling between the two junctions, 
i.e. mutual capacitance between the normal wires, prop- 
agating modes in the superconducting wire^^ and similar 
effects. 

Let us first " turn off" electron-electron interactions by 
taking the limit l/Gf^2 ~^ and, hence, setting (pr = 
eVrt. Defining the non-local noise spectrum as Sri{uj) — 
2j dtS^j' exp{~iuj{t - t')), from Eq. ^ we obtain 

5i2H = -2Gi2(2-/?i-/32)T^(w,0) (35) 

- 2Gi2PiW{lu, 2Fi) - 2Gi2/32VF(w, 2F2) 

+ Gi2j+W{w, Vi + V2) - Gi2l-W{uj, Vi - F2), 



IV. INTERACTION CORRECTION TO THE 
CURRENT 

Now we again "turn on" electron-electron interactions 
and evaluate the current /i across the first barrier. Solv- 
ing Eqs. (|32|) perturbatively in l/g^ <C 1, in the lowest 
non-trivial order in this parameter we get 



h = GiiVi - G12V2 - (6) 



(39) 



Here the average (fi) does not vanish since according to 
Eqs. (I?71) . (1^51) the noise ^1 depends on the phases V'1,2, 
which, in turn depend on ^1.2 by virtue of Eqs. (j32p . 
Hence, we obtain 



where 



W(,u:,V)^\Y.^u±eV)coth'^^^. (36) 



(Ci) = {5ipi d^i/dipi) + {Sip2 d^i/dip2 



(40) 



where the phase fluctuations Sifj. are found from Eqs. 
(l32l) and read 



This is the first key result of the present paper. Eqs. (|35|) , 
([M)l fully determine non-local shot noise correlations in 
NSN structures at subgap energies T,uj,eVr <^ |A|. In 
the tunneling hmit T12 < 1 we have 7± = 4:/TiT2 > 
/3i^2 — 1 and Eq. (|35|) is dominated by the last two 
terms which further reduce to the result^ in the limit 
w — )■ 0. One of these contributions ex 7_ is due to EC and 
describes negative shot noise correlations while the other 
one oc 7+ comes from CAR and accounts for positive 
cross-correlations of fluctuating currents. Provided one 
of the voltages, say Vi, equals to zero, these EC and CAR 
terms exactly cancel each other for any V2, i.e. 5i2(0) 
tends to zero in the tunneling limit. 

At higher transmissions ti_2 the value 7_ becomes neg- 
ative thus implying positive cross-correlations produced 
by direct electron transfer at such values of T1.2 in ad- 
dition to CAR. In the limit of fully transparent barriers 
ri_2 = 1 one has /3i,2 = 7_ = 0. Then Eq. (|35]) yields 



5i2(0) = -8rGi2 



- 2eGu{Vi + V2) coth ^^^'^^^'^ (37) 



At T — ^ only positive cross-correlations due to CAR 
survive whereas no direct electron transfer contribution 
to shot noise occurs for fully open barriers. Accordingly, 
no non-local shot noise is expected in this case at Vi — 
— V2 and T — > 0. We also note that noise correlations 
in clean NSN systems and in the specific limit Ti = T2, 
V\ = V2t T — uj = Q have been recently studied in Ref. 
I25I . The result of this paper is consistent with our Eq. 
([55)) in the corresponding limit. 

For completeness we also provide the expression for the 
local noise 5*11(0;), which is given by the Fourier trans- 
form of Eq. ([27]) and reads 

5ii(w) = 2Gii [{\ ~ h)W{uj,Q) + liiW(uj,2Vi)] 

+ 4Gi2 [aiW{uj, 0) - TiiW{uj, 214)] 

+ Gi2i4W{oj, Vi + V2) + Gi2KlW{ijj, Vi - y2)(38) 



5ipr{t) 



1 _ p-{i-t')/'rR.C 

dt' 7^ Ut')- (41) 



Here we have assumed G12 <C Grr ^ G'^^ and intro- 
duced the i?G— time trc — Cr/Gf^. Substituting this 
expression into Eq. PO]) we find 



(6) 
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dipr{t) 



(.42) 



ip^=eVrt 



Making use of Eqs. p7|. (|28p and performing the time 
integral in Eq. (|42p we get the current through the first 
junction in the form 



h — Gii^i — G12V2 
2Gii/3i - 4Gi2r?i 



Fo{2Vi) + ^^^^Fo{2V2) 



.91 92 

-SG+Fo{Vi + V2) - 5G-Fo{Vi - V2), 

where dG± — G12 (k^ /gi + j±/g2) and 
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.(44) 



Accordingly the local and non-local differential conduc- 
tances read 



dh ^ 4Gii/3i-8Gi277i ^,,^,, 
- 5G+F{Vi + V2) - SG^F{Vi - V2), 



(45) 



and 



dV2 - ~^'' 



1 _ ^Fi2V2) 

92 



■5G+F{Vi + V2)+5G-F{Vi-V2), (46) 
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FIG. 2: Local (a,c) and non-local (b,d) differential conduc- 
tances defined respectively in Eqs. 1)45^ and ((46]). The param- 
eters of the system are: T = 20 mK, Gn = 1 mS, G12 = 10 
fiS, gi ~ 92 = 516. The dependencies (a) and (b) are plotted 
at small transmissions T\ = 0.063, T2 = 0.11, while the graphs 
(c) and (d) correspond to fully open barriers Ti = T2 = 1. 
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and '^{x) is the digamma function. Eqs. (H5|) . P5|) to- 
gether with Eq. (|T7)) define the conductance matrix of 
our NSN device and represent the second key result of 
our paper. 

In the interaction correction to the local conductance 
in (|45p we recover the Coulomb blockade term*' ex /3i 
and, in addition, three non-local contributions. The first 
of them ex 771 enhances the conductance, while the second 
one ex 5G-i^ provides additional Coulomb suppression of 
dli/dVi. The last term ex 5G- can be both positive (at 
Ti,2 ^ 1) and negative (at bigger ri,2) implying the ten- 
dency to Coulomb anti-blockade in the latter case. The 
first term ex P2 in Eq. (|^^ has an opposite sign as com- 
pared to G12 (thus implying Coulomb blockade), while 
the second one ex 5G+ yields Coulomb anti-blockade. Fi- 
nally, the third ex 5G- tends to suppress or enhance the 
absolute value of the non-local conductance respectively 
for 5G- > and SG- < 0. 

The origin of each of the terms in both interaction 
corrections can easily be identified from the correspond- 
ing shot noise correlators (P7|) . (PSJ) and ([55)1 illustrat- 
ing again a fundamental relation between shot noise and 
Coulomb effects in electron transport. This relation turns 
out to be considerably more complicated than in the lo- 
cal case. In the tunneling limit Ti,2 ^ 1 in Eq. ps)) 



the non-local terms add up to the local one and dli/dVi 
evolves from a typical Coulomb blockade V-like depen- 
dence at small V2 towards a new W-like one (with ex- 
tra minima at Vi — ±¥2) at higher V2 (Fig. [5^). In 
Eq. (|46|) . in contrast, the last two terms exactly cancel 
each other for Vi — >■ and any V2 since 6G+ ~ SG-. 
This cancellation has the same origin as that of EC and 
CAR contributions to shot noise discussed above. For 
nonzero Vi and V2 the last two terms in Eq. (|46)) do not 
cancel anymore and the curve dIi/dV2 approaches the 
S-like shape with maximum at Vi = V2 and minimum at 
Vi — —V2 (Fig. ^jp). In this case the interaction term 
ex SG± ex T1T2 {l/gi + I/32) can exceed G12 ex T^T^ 
and, hence, dIi/dV2 can change its sign. For fully open 
contacts with Ti_2 = 1 we get /3i.2 = and 5G- = 0, i.e. 
only CAR terms containing SG+ = 2Gi2 (I/51 + l/<?2) 
survive in Eqs. (|45|) and (|46|) implying Coulomb blockade 
for local conductance dli/dVi and anti-blockade for non- 
local conductance \dIi/dV2\ in this limit (Figs. [2t,d). 

Finally, we would like to note that in some cases non- 
linearities in both local and non-local differential con- 
ductances caused by Coulomb interaction may be com- 
bined with the zero bias anomalies resulting from the 
proximity-enhanced electron interference in diffusive nor- 
mal leadsiii22r— . In this paper we disregarded this effect 
for the sake of simplicity. In practice it implies that here 
we considered the system with weakly disordered or suf- 
ficiently thick normal leads and sufficiently resistive bar- 
riers. If needed, zero-bias anomaly effectsi^i^i"— can be 
included into our analysis in a straightforward manner. 



V. SUMMARY 

In this paper we developed a theory elucidating a non- 
trivial physical relation between shot noise and Coulomb 
effects in non-local electron transport in NSN structures. 

We evaluated non-local current-current correlators in 
such systems at arbitrary transmissions of NS interfaces 
and arbitrary frequencies, Eqs. (P5|) . (P5|) . This re- 
sult demonstrates that positive cross-correlations in shot 
noise increase with increasing interface transmissions and 
dominate the result for fully open barriers in which case 
only CAR contribution survives. Positive noise cross- 
correlations in NSN structures have been convincingly 
demonstrated in recent experiments^^, while no or weak 
negative cross-correlations have been observed. This pic- 
ture would qualitatively correspond to the case of highly 
transmitting interfaces, cf. Eq. ([57)) . Note, however, that 
interface transmissions in experiments^S are reported to 
be rather small, in which case one would expect EC- 
induced negative cross-correlations to dominate the re- 
suh at Vi = -^2- 

Turning to the effect of electron-electron interactions 
on non-local electron transport we would like to empha- 
size several important new features demonstrated within 
our analysis. One of them is that in the tunneling limit 
almost no effect of Coulomb interaction on non-local con- 



ductance is expected if one of the applied voltages, Vi 
or V2, equals to zero. This effect is directly related to 
the cancellation between EC and CAR contributions to 
shot noise in the corresponding limit. For nonzero Vi 
and V2 no such cancellation exists anymore and the non- 
local conductance dIi/dV2 approaches the S-like shape 
being enhanced at Vi « V2 and partially suppressed at 
Vi w — V2, see Fig. [5]3. Both these features have a 
clear physical interpretation. Indeed, at Vi ~ —V2 nega- 
tive cross-correlations due to EC dominate non-local shot 
noise leading to Coulomb blockade of non-local conduc- 
tance while at Vi ~ V2 positive cross-correlations due 
to CAR prevail and Coulomb anti-blockade of non-local 
transport is observed. At higher interface transmissions 
only Coulomb anti-blockade of non-local conductance re- 
mains (Fig. [2Ji), which is again related to CAR- induced 
positive cross-correlations in shot noise. 

It is interesting to point out that S-like shaped 
non-local signal predicted here was indeed observed in 
experimentsiii^. A good agreement between our theory 
and the results^^ argues in favor of electron-electron in- 
teractions as a physical reason for the observed feature. 
Some of the features similar to those predicted here have 
also been observed in experiments^^. More experiments 
on both non-local shot noise and non-local electron trans- 
port would be desirable in order to quantitatively verify 
our predictions. 
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Appendix A: Parameters of the action 

For the sake of completeness, let us present general 
expressions for non-local conductance G12 and the pa- 
rameters K^ and 7± entering in the effective action (j27l) , 
(1281). We have 



G12 - — tr [{X-i - Xi)KR{Xi - X2)Ka] , (Al) 
Zn 



(e^FiFs/TrGia) tr [2KRXiKAXi 

+ 2TlKRX,KAiX, ~ X2)KrX^Ka{Xi - X2) 

- 2iFi(Xi - X2)KrX^KrX^Ka 

+ 2iTiXiKa{X^ - X2)KrX^Ka] , (A2) 

(e^FiFa/TTGia) tr [RrX^RrX^ + KaX^KaXi 

- T{r2KRXiKA{Xi - X2)KRXiKA{X^ ~ X3) 

- Y^V2KrX^Ka{X^ - X3)KRXiKA{Xi - X2) 

+ ir2iXi - X3)KrXiKrX^Ka 

+ iTi{Xi - X2)KrX4KrXiKa 

- iT2XiKA{Xi - X3)KrX^Ka 

- iTiXiKA{Xi - X2)KrXiKa] ■ (A3) 



7+ 



The parameters k^ and 7_ are defined by Eqs. (jA2| and 
(jA3P with interchanged matrices X^ <->■ X^ . The remain- 
ing parameters a^ and 77^ in Eqs. (P?]). (|28l) are defined 
in a similar manner. These parameters are less impor- 
tant for our consideration and we omit the corresponding 
expressions for the sake of brevity. 
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